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Zusammenfassung 

Notwendiges und hinreichendes Kriterium dafür, dass eine Fixpunktiteration vom Typ 

|6R" X "; x m ,h e R"; m £ No, n G N 

unabhängig vom Anfangswert gegen einen Fixpunkt x* £ R" konvergiert (für den dann: 
x* — $ x* + h erfüllt ist) ist bekanntlich, dass für alle Eigenwerte von J? : Ai, A; (l < n) 
gilt: | Xi |< 1. i G {1,...,/} 

In dieser Arbeit wird ein neues Verfahren angegeben, das es gestattet die Konvergenz der 
Fixpunktiteration auch für den Fall zu erzwingen, dass die Eigenwerte von $ (der " Iterati- 
onsmatrix") betragsmäßig größer als eins sind. 

Im ersten Kapitel dieser Arbeit wird ein neues Konzept vorgestellt, das es erlaubt, eine vor- 
gelegte Fixpunktiteration so zu erweitern, dass jedes Eigenwerttupel der neu entstehenden 
Iterationsmatrix eingestellt werden kann. Das Eigenwerttupel der ursprünglichen Iterations- 
matrix $ muss aber noch bekannt sein. Ein notwendiges und hinreichendes Kriterium für die 
Durchführbarkeit dieses Konzepts wird angegeben. Aus den errechneten Ergebnissen kann 
sehr einfach auf die Lösung des Ausgangsproblems zurückrechnet werden. 

Im zweiten Kapitel wird dieses Konzept dann so erweitert, dass man auf die Kenntnis des Ei- 
genwerttupels von J? verzichten kann: Die Verwendung von Optimierungsverfahren für lineare 
dynamische Systeme erzwingt eine Iterationsvorschrift, bei welcher der spektrale Radius der 
Iterationsmatrix kleiner eins (aber ansonsten nicht genauer bekannt) ist. 

Das dritte Kapitel liefert schließlich eine Iterationsvorschrift, bei der der spektrale Radius 
vorgeschrieben werden kann, ohne dass das Spektrum der Iterationmatrix $ bekannt sein 
muß. Insgesamt wird sich zeigen, dass sich damit das Einsatzgebiet von iterativen Rechenver- 
fahren zur Lösung von linearen Gleichungssystemen bedeutend erweitern läßt. 
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1 Ein neues Konzept für die iterative Berechnung der Lösung 
eines linearen Gleichungssystems 

Für die Überführung eines linearen Glcichungssystems vom Typ 



A ■ x 



in die Matrixdifferenzengleichung 



m e No, x — R" 



wird hier beispielhaft das Jacobische Iterationsverfahren verwendet. Aber auch für andere Iterati- 
onsverfahren (wie z.B. das Gauß-Seidelsche Verfahren) bleibt das hier vorzustellende Lösungskon- 
zept gültig. 



1.1 Eine Erweiterung des Jacobischen Iterationsverfahrens 

Ausgehend von 

A-x = b ; A G R" x ™ b,xGK n 



(1) 



erhält man bekanntlich nach Einführung von: 

Ll = {hj)i,j=l,...,n h,j 
H= (Uij)i,j=l,...,n 
— = {hj)i,j=l,...,n 

für (1) die Darstellung: 





% > j 


: 


sonst 


a i,3 


: i <j 





: sonst 




: i = j 





: sonst 



d i,j = 

(L + D + U)x = b 



(2) 



und weiter (wenn JD invertierbar ist) 



x= -IT 1 ■ (L + U) -x + IT 1 -b (3) 

Kürzt man ab: 

$ = -Dr 1 ■(L + U) (4) 

h = D~ 1 -b (5) 

so kommt man auf die zugehörige Fixpunktform 

x = $-x + h (6) 

Davon ausgehend wird die Iterationsvorschrift (Matrixdiffcrcnzcnglcichung) 

x m+1 = ± ■ x m + h; meN ,x = R™ (7) 

formuliert. Man definiert: 

Definition 1.1 Das Iterationsverfahren (7) ist konsistent zur Matrix A, wenn für alle b G R™ 
der Term A^ 1 ■ b ein Fixpunkt für dieses Iterationsverfahren ist. 
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Bekanntlich ist die Itcrationsvorschrift (7) konsistent zur Matrix A, wenn für den Spcktralradius 
von $ gilt [Mci] : 

P®<1 (8) 
Für den Fall p ($) > 1 werde die Iterationsvorschrift (7) jetzt wie folgt erweitert: 

x m+1 = (ß_ — h- k T ^j -x m + h; Xq = x ] k <E R" (9) 

=* 

Anmerkung 1.1 Die Differenzengleichung (9) wird in dieser Arbeit das erweiterte Iterationsver- 
fahren genannt. 

Durch die Forderung nach einem Eigenwerttupel für $ soll der Vektor k festgelegt werden. 
Diese Forderung beschreibt das (in der Regelungstechnik seit langem verwendete) Konzept der 
Zustandsregclung. Unter welchen Bedingungen diese Forderung erfüllt werden kann, besagt der 
nachfolgende Satz [Hartl]: 

Satz 1.1 Notwendig und hinreichend dafür, dass für eine Matrix l> = <& — h ■ fc T durch geeignete 
Wahl von k jedes beliebige Eigenwerttupel eingestellt werden kann ist, dass: 

rgiä^-h,---^! 1 ' 1 -h) =n 

erfüllt ist. 

Rechenverfahren zur Berechnung eines geeigneten k - Vektors nachdem 
die Eigenwerte Ai , • • ■ , A„ der Matrix $ bekannt 

und 

die Eigenwerte Ai, • • • , A„ der Matrix J> gefordert 

werden, sind seit langem bekannt. Sie laufen auf einen Koeffizientenvergleich der zughörigen cha- 
rakteristischen Polynome (von $ und $ ) und einen Basiswechscl im R" hinaus. Um aber diejenige 
Transformationsmatrix zu berechnen, die diesen Basiswechsel vermittelt, ist die Lösung eines wei- 
teren n x n - Glcichungssystcms erforderlich. [Focll] 



Damit erscheint das gesamte, neue Konzept für die iterative Berechnung der Lösung eines linearen 
Glcichungssystcms hinfällig. Der Aufwand ist normalerweise zu groß: damit das Glcichungssystem 
(1) in der Darstellung (9) iterativ gelöst werden kann, müßte ja vorher ein anderes Glcichungssy- 
stem von derselben Ordnung gelöst werden. Hinzu kommt die Berechnung des Eigcnwerttupels von 
$. Immerhin sollen einige bemerkenswerte Aspekte für dieses neue Konzept vorgestellt werden. 
Damit rechtfertigen sich weitere Modifikationen die schließlich auf ein leistungsfähiges Lösungs- 
konzept (vgl. Folgekapitel) führen. 

Von den neuen Aspekten seien genannt (diese gelten, wenn die Bedingung von Satz 1.1 erfüllt ist): 

1. Unabhängig vom Spcktralradius der Matrix $ kann man eine gegen den Fixpunkt von (9) 
konvergente Iteration erzeugen. 

2. Indem man von $ verlangt: Ai = 0, • • • , A„ = 0, konvergiert die Iteration (9) nach genau 
"n" - Iterationsschritten und das Ergebnis ist exakt. 

3. Die Rückrcchnung von dem unter (9) errechneten Fixpunkt in den Fixpunkt von (7) und 
damit in die Lösung von (1) ist elementar. 

Bevor diese neuen Aspekte anhand eines kleinen Rechenbeispiels demonstriert werden, soll der 
unter Punkt 2 genannte Aspekt hervorgehoben werden. 
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1.2 Minimale Anzahl von Iterationen 

Die Bedingung von Satz 1.1 sei erfüllt, k kann dann für jedes geforderte Eigenwerttupel von $ 
(also von Ai, • • • , A„) berechnet werden. Insbesondere auch für Ai = 0, • • • , A„ = 0. $ ist dann 
nilpotent. Dieser Fall soll jetzt für die $ - Matrix der Differenzengleichung (9) vorliegen: 

Der zugehörige Fixpunkt sei x: 

x = Ö_-x + h 

Subtrahiert man die letzten beiden Gleichungen voneinander, so erhält man 

iE m _l_i — x_ = $ • (x_ m — x) ; x_ g — x = x_ q — x (10) 

und mit 

Sm • £m — 3L (H) 

endlich 

Sim+l = ^ ' £-m > ^0 = 3k) ~ 3L (12) 
Rückwärtseinsetzen der einzelnen Gleichungen aus (12) liefert 

Si = t • e 
£2 = $ • e 



Wegen der vorausgesetzten Nilpotenz von $ ist e n = und folglich i n = i. Nach m = n Schritten 
wird der Fixpunkt also exakt erreicht. 

Im anschließenden Beispiel wird noch gezeigt, wie sich aus diesem Fixpunkt x für das erweiterte 
Iterations verfahren die Lösung x für das Ausgangsproblcm (1) berechnet. 

1.3 Beispiele 

1.3.1 Sämtliche Eigenwerte von J> liegen bei Null 

Für das entsprechend (1) zu lösende, lineare Gleichungssystem lauten Matrix und Vektor 

'1 N 

A = I 29 42 511 und b = I 266 I , mit Lösungsvektor x = I 2 



59 


-63 








29 


42 


2 


und 


b 


30 


31 


-67/ 








Die zugehörigen Matrizen entsprechend (7) sind 



0,000.000.000 1,067.796.610 0, 033.898.305 \ / -1, 237.288. 135.593N 

$=|-0,690.476.190 0,000.000.000 -1,214.285.714 und h = 6,333.333.333.333 
0,537.313.433 0,462.686.567 0,000.000.000 / \ 1,537.313.432.836, 

Da Ai = 0, 235.428 + i ■ 1, 202.990, A 2 = 0, 235.428 - i ■ 1, 202.990, A 3 = -0, 470.856 die Eigenwerte 
von $ (und damit |Ai|, | A.2 1 > 1) sind, wird das Iterationsverfahren (7) nicht konvergieren. 
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Verlangt man nunmehr, dass $ (vgl. (9)) einen Dreifacheigcnwert bei Äi = 0, A2 = 0, A3 = 
habe, dann ergibt sich (nach Kocmzicntcnvcrglcich der charakteristischen Polynome von $ und $ 
und anschließendem Basiswechsel im R" [Focll]) ein Vektor 

k T = (-0,117.899 0,025.369 -0,199.404) 

Nach drei Schritten erreicht die Itcration (9) den Vektor 



±3 



2,988.423> 

5,976.846 

8,956.270; 



(Der Startvektor war in diesem Fall der Nullvektor.) Zufolge 1.2 ist (unabhängig vom Startvektor) 
x 3 der Fixpunkt x der Iteration (9). 

Den originalen Lösungsvektor "x" erhält man, wenn man, ausgehend von (9), den Wert für x 
einsetzt 

£ = (t - h ■ k J 

auflöst (vgl. (4) und (5)) 

x = $ • x — h-k T -x + h 
= -DT 1 -(L + U)-x~ LT 1 ■ b ■ k T ■ x + LT 1 ■ b 
mit D auf beiden Seiten der Gleichung durchmultipliziert 

D-x = -(L + U) ■ x-b- k T ■ x + b 

umstellt und rücksubstituiert: 

A-x = b — b-k T -x 

Dann erhält man 

x = A- 1 -b-A- 1 -b-k T -x 

( T 

und endlich 

i (13) 



1 — k T ■ x 



Berechnet man schließlich mit diesem Lösungsvektor den Restfehler r := A ■ x — b, so erhält man 
für das vorgelegte Beispiel Restfehlerkomponenten von der Größenordnung 10~ 14 . 



1.3.2 Die Eigenwerte von $ liegen bei Ai = 0, 1 , A 2 = 0, 2 und A3 = 0, 3 

Hierfür lautet 

k T = (-0,119.694 -0,072.227 -0.189.066) 

Um die Genauigkeit von 1.3.1 zu erreichen sind 45 Iterationsschrittc notwendig. Es ist 

/ 5, 929.411 \ 
i 45 = 11,8 58.823 
\ 17,788.235/ 

Der Startvektor für dieses Iterationsergebnis ist xq = (1 1 1) T . Nach Rücktransformation (13) 
ergibt sich der Lösungsvektor mit Restfehlerkomponenten wieder von der Größenordnung 10~ 14 . 
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2 Anwendung von Optimierungsverfahren 



Bestimmte Optimierungsverfahren für lineare Regelsysteme liefern ebenfalls einen Vektor "fc" der 
einen Spektralradius p (jh^J < 1 erzwingt. Es ist naheliegend, diesen Vektor "fc" auch für das er- 
weiterte Itcrationsverfahren (9) einzusetzen. Diese Optimierungsverfahren benötigen keine Berech- 
nung von Matrixcigenwerten. Die Grundlagen zu dem hier vorzustellenden Optimicrungskonzcpt 
sind z.B. in [Ludyk] dargestellt. Eine typische Aufgabenstellung und die zugehörige Lösung für 
die Optimierung linearer, dynamischer Systeme wird jetzt vorgestellt. Anschließend wird der, für 
das Optimierungsverfahren gefundene, Vektor "fc" in unserem Iterationsverfahren verwendet. 

2.1 Das Optimierungsproblem 

Gegeben sei für $ G R nx ™, h G R™, {u m } mG N mit u m G R 
die Iterationsvorschrift für die Folge {x m } m eTSio m it Um £ 

x m +i = $ ■x m + h-u m (14) 

Weiterhin seien gegeben: 

• die positiv defmite Matrix S_ G R nxrl 

• die Zahl r > 

• ein Vektor c G R™ mit der Eigenschaft, dass gilt: 



Rang 



= n und C'C T = Qe R™*" 1 ist positiv semidefinit (15) 



Gesucht ist für beliebig vorgegebens N G N eine Vektorfolge {k m } mS {o,i,2 n-i}, sodass 

nach Bildung von 



z.m+1 ~ ^ ' —m + ^ ' u m\ z.m=o — 2-0 (16) 

u m = — hm ' Z.m> (l'O 

das quadratische Gütekritcrium 

2 



1 1 N ~ 1 

^ = x^-5-z Ar + -^z^(c.c T )z m + r-^ JVeN (18) 



m=0 

minimal wird. 



Der Wechsel der Variablenbezeichnung von "x m " nach "z m " wird notwendig, da nach Einsetzen 
von (17) in (16) im Gegensatz zu (9) nunmehr eine homogene Matrix - Diffcrenzengleichung 
entsteht. Man beachte weiterhin, dass in (16) mit variablem "k m " gerechnet wird. 



Die Lösung dieser Optimierungsaufgabe lautet unter der Voraussetzung, dass zusätzlich Satz 1.1 
gilt [Ludyk]: 
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Satz 2.1 Für jedes N G N gilt: Für das mathematische Modell (16)-(18) in Verbindung mit der 
Forderung: 

Jn = Min 

berechnet sich der "optimale Vektor" zu 



T 



1 



k 

—in . ,T D 

r + h ■ P m+1 ■ h 



(h T ■ P m+1 ■ t); me{0,l,2,..,iV-l} (19) 



Hierbei ist P_ m £ R, nx ™ efjg Lösung der Matrix - Riccati - Differenzengleichung 



P m = Q + tL 



T 



T 1 

p _ p hh P • 

i_m+l £-m+l'-L!± £-m+l ,T 



h E m +ih 



#; me{l,2,...,iV-l} (20) 



mii dem Endwert P_ N = S_. 



Die Matrix - Riccati - Diffcrenzengleichung ist also "rückwärts" zu lösen. Aus der so entstehenden 
Folge {P_n,P.n-i> ' • • errechnet sich {k N _ 1 ,k N _ 2 , ■ ■ ■ ,k }. 

Der sich für " J^r" ergebende Wert ist für unsere Betrachtungen nicht weiter interessant. 
Wichtig hingegen ist, daß die damit verbundene Rckursion (die homogene Matrix - Diffcrenzen- 
gleichung) 

unabhängig vom Anfangsvektor z gegen den Grenzwert limjv^oo z N = konvergiert. (Das ist 
zumindest plausibel, da wegen der vorausgesetzten positiven Dcfinithcit von S_ und der Scmidcfi- 
nitheit von Q ein Minimum für " Jn" nur angenommen werden kann, wenn liniTv-yoo z N = gilt. 
Man beachte, daß für "N" keine weiteren Annahmen als TV G N gemacht wurde und dass folglich 
"N" beliebig groß sein darf. Tatsächlich kann man zeigen, dass aufgrund der vorausgesetzten Ei- 
genschaften von Vektor c, vgl. (15), die positive Semidefmitheit von Q hierfür ausreicht. [Ludyk]) 



Anmerkung 2.1 In der nach (19) berechneten Vektorfolge {fc m } mg | 1 2 n-i} konvergiert für 
N — > oo und unter der Voraussetzung (15) der Vektor k gegen einen Grenzwert. Dieser werde 
wieder mit "k" bezeichnet. Der Grenzwert k ergibt sich auch als Lösung der diskreten algebraischen 
Matrix - Riccati - Gleichung: 



P = Q + * 



P-P hh 1 P. 



1 



r + h T P h 



(21) 



zu 



e = 



i 



h 2 P h 



h 1 P$ 



(22) 



Obwohl für die Lösung der diskreten algebraischen Matrix - Riccati - Gleichung schon Lösungspro- 
zeduren bereitgestellt sind (z.B. der Befehl "DARE" unter Matlab), empfiehlt es sich im Hinblick 
auf das betrachtete Problem (iteratives Lösen von linearen Glcichungssystcmcn) , doch besser mit 
der Matrix - Riccati - Diffcrenzengleichung (20) zu rechnen: alle uns bekannten Lösungsprozedu- 
ren für die Lösung der diskreten algebraischen Matrix - Riccati - Gleichung verwenden intern die 
Berechnung von Matrixcigcnwcrten, was bei Glcichungssystcmcn hoher Ordnung wieder auf Nu- 
merikproblcmc führen kann. Durch Einsetzen von P_ m in (21) kann man die Qualität der iterierten 
Lösung überprüfen. Sei: 



A = Q + J? 



P, 



P„ h h 1 P„ 



1 



h T P m h 



und liege 1 1 A 1 1 unterhalb einer hinreichend kleinen Schranke, so kann P_ m durch P_ ersetzt werden. 
Damit berechne man fc T entsprechend (22). 
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Anmerkung 2.2 Mit diesem k - Vektor konvergiert aber auch die Iteration (9): 
2Lm+l = (ß?- ~ h. ' k ) • 2L m + hü 3±m=V> = SLq 

gegen einen Fixpunkt. 

Bevor weitere Anmerkungen zu den Matrizen S_, Q und P_ m sowie zu "r" gemacht werden, sollen 
die bisherigen Ergebnisse an einem Beispiel demonstriert werden. Wir greifen wieder die A und 
die $ - Matrix sowie den b und h - Vektor von Beispiel (1.3) auf. 

2.2 Weiterführung des Beispiels 

Mit den recht willkürlichen Werten für 



„T 



c 



(14 5) 



und folglich 



l 4 5 
0=4 16 20 
\5 20 25 y 

sowie mit 

S_ = E - Einheitsmatrix und r = 0, 5 

erhält man für N = 100 Iterationen, beginnend mit Pioo = Einheitsmatrix mit Hilfe der 
Gleichung (20) 



'1,482.773.652 4,200.001.512 5, 222.175.153 \ 
4,200.001.512 17,345.911.012 19,930.235.464 
,5,222.175.153 19,930.235.464 25,132.197.634/ 



'1,035.281.548 4,023.046.495 5, 014.945.202 ' 
P 2 «( 4, 023.046.495 16,060.198.199 19,996.629.128 
,5,014.945.202 19,996.629.128 25, 018.325. 346> 



'1,035.281.548 4,023.046.495 5,014.945.201' 
Pi « | 4, 023.046.495 16,060.198.200 19,996.629.128 
,5,014.945.201 19,996.629.128 25, 018.325. 346 1 



Für die Vektoren k^ ergibt das mit der Gleichung (19): 



km ~ (-°) 079.698.253 - 0, 013.703.479 - 0, 173.741.413) 
km ~ 007.838.204 0, 102.143.449 - 0, 153.895.906) 



k{ » (-0, 002.607.656 0,106.247.305 -0,151.791.016) 



fco « (-0, 002.607.656 0, 106.247.305 - 0, 151.791.016) 



8 



Und kompakter (jetzt mit gerundeten Zahlenwerten): 



m 


p 


h 


100 


1,000 0,000 0,000 0,000 1,000 0,000 0,000 0,000 1,000 




99 


1,483 4,200 5,222 4,200 17,346 19,930 5,222 19,930 25,132 


-0,079.698 -0,013.703 -0,173.741 


98 


1,672 4,629 5,258 4,629 17,168 20,160 5,258 20,160 25,119 


-0,007.838 0,102.143 -0,153.896 








1 


1,035 4,032 5,015 4,023 16,060 19,997 5,015 19,997 25,018 


-0,002.608 0,106.247 -0,151.791 







-0,002.608 0,106.247 -0,151.791 



Die nachstehende Abbildung 1 zeigt den Verlauf der Komponenten des iterierten fc m -Vektors: 

k1 = o,k2 = *,k3 = + 

0.15 



0.1 



0.05 



1— 

o 



-0.1 



-0.15 



-0.2 




1 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I | 1+4 



10 20 30 40 50 60 70 

Iterationsnummer (Rückwärtsiteration) 



80 



90 



100 



Abbildung 1: Darstellung der Komponenten des k m - Vektors. Die Konvergenz gegen einen Grenz- 
wert ist unter Vorlage der Voraussetzungen von Satz 2.1 und von (15) gesichert. 



ij ist eine Näherungslösung von P_ entsprechend (21): 

/l, 035.281.547.600 4,023.046.494.700 5,014.945.201.200' 

P = 4,023.046.494.700 16,060.198.198.000 19,996.629.128.000 

V 5, 014.945.201.200 19,996.629.128.000 25,018.325.346.000; 
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und k Q eine Näherungslösung von (22): 

k T = (-0, 002.607.655.915 0, 106.247.305.105 - 0, 151.791.015.721) 

Mit diesem k- Vektor und dem Startvektor x£ = (1, 1, 1) liefert die Iteration (9): 

0, 758.352.051. 190\ 
1,257.479.544.079 
3,074.126.036.400; 



0,876.097.560.719> 

2,157.717.086.365 

3,003.581.738.766; 



1, 


140. 


.376. 


318 


225\ 


1, 


890, 


942. 


195 


607 


4, 


622. 


,734. 


949 


,616/ 



1. 


136 


.239 


.118 


,708\ 


'2: 


798 


.412 


.723. 


144 


3, 


895 


.441 


.810. 


,180/ 



±250 



1, 


325 


.309 


.503. 


,590^ 




2, 


650 


.663 


.291 


,458 


| 2Lioo 


3, 


976 


.119 


.363 


.656) 




1. 


325 


.356 


.617. 


75D, 




2. 


650 


.713 


.325 


.513 


| ^250 


3, 


976 


.069 


.853, 


,246) 





0, 999.954.870. 724^ 

1,999.943.154.231 

3,000.008.611.948,1 



1, 000.000.000. 002\ 

2,000.000.000.011 

2,999.999.999.999/ 



(Bei der Berechnung dieser x- Vektoren, wenn statt "fc" der Vektor "fc " verwendet wird, hat man 
für x 10 und für x 25 eine Ubereinstimmung in den ersten 9 Nachkommastcllen; bei x wo und x 250 
jedoch volle Ubereinstimmung.) 

Hierbei beträgt der Rcstfehler r m := A ■ x m — b: 



-100 



0,901.432 10~ 3 \ / -0,619 10" 

-3,257.061 10" 3 r 25 o= 0,486 10~ 9 I (23) 



3,963.873 lO" 3 / V 0,482 10" 9 



Für kleine m- Werte sind die Restfehler r m aber noch " sehr groß" . Die langsame Konvergenz der 
Folgen (i m ) mgNo und (^ m ) m6No wird durch die Eigenwerte von J> bedingt. Diese müssen ja nicht 
mehr vorgegeben werden - die Lösung des Optimierungsproblems garantiert lediglich |Aj| < 1. 
Tatsächlich liegen diese mit 

|Ai| = 0,900.828.330.426, |Ä 2 | = 0,900.828.330.426, |Ä 3 | = 0,000.430.484.178 

zum Teil nahe bei eins. Die Problematik der gezielten Beeinflussung des Eigenwerttupels (ohne 
dessen explizite Berechnung) wird in Kapitel 3 aufgegriffen. 
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2.3 Abschließende Bemerkungen zum Optimierungsverfahren 



Für das in 2.1 vorgestellte Optimierungproblem wird eine Lösung angegeben, die jetzt auch die 
Berechnung eines linearen Gleichungssystems (1) durch zwei Iterationsverfahren erlaubt: 

1. Iterative Berechnung vonPj durch Rückwärtslösen der Matrix- Riccati-Diffcrcnzengleichung. 
Anschließend Berechnung von k . 

2. Mit k = k das erweiterte Fixpunktverfahren (9) durchführen. Anschließend Berechnung des 
Lösungsvektors "x" durch Rücktransformation entsprechend (13). 

Da jedoch von den Eigenwerten von J> nicht mehr gesichert ist, als: | Ai |< 1 i € {1, / < n}, 
so kann noch keine Genauigkeitsangabe (a-priori-Fehler und a-postcriori-Fchler) für das erweiterte 
Iterationsverfahren gemacht werden. Eine solche Angabe gelingt aber, da mit Hilfe des vorgestell- 
ten Optimierungsverfahrens auch der betragsgrößte Eigenwert von J> abgeschätzt werden kann. 
Dies wird im nächsten Kapitel ausgeführt. 



3 Ergänzungen zum quadratischen Gütekriterium 

Die Auswahl der Matrizen "S", des Vektors "c" und der Zahl "r" scheinen zunächst im Hinblick 
auf die vorgelegte Aufgabenstellung: "iterative Lösung eines linearen Gleichungssystems" viele 
Freiheiten zu lassen. Die Wahl von "S" erscheint unkritisch, da die Matrix - Riccati - Differen- 
zengleichung unabhängig von "S" gegen die Lösung der diskreten algebraischen Matrix - Riccati - 
Gleichung konvergiert. 

Zusätzlich zu der Tatsache, dass für den Spektralradius von 

erlauben die Wahl von "c" und "r" eine genaue Platzierung der Eigenwerte von $ im Innern des 
Einheitskreises der komplexen Zahlenebene [Hart2] . Da dazu aber auch die Eigenwerte von $ zuvor 
bekannt sein müssen, wird (im Hinblick darauf, dass das hier vorzustellende Iterationsverfahren 
für Gleichungssysteme hoher Ordnung angewendet werden soll) dieser Aspekt nicht weiter verfolgt. 

Immerhin kann der Spektralradius p I abgeschätzt werden. 



3.1 Schranken für den Spektralradius p(j$ij. 

Ausgehend von Iteration (14) mit $ —> ^ • $ =: $ w und h^^-h~:h w (w > 0) 
und 

2Lm + l = 3Lw ' 3Lm + thti ' u m i 2Lm=0 = £o (24) 

werde eine Folge {u m } m gN gesucht, sodass mit u m — —k^ ■ x m das quadratische Gütekriterium 
(18) minimal wird. Sei fc* der "optimale Vektor" der sich als limAr^oofcQ oder als Lösung der 
diskreten algebraischen Matrix - Riccati - Gleichung ergibt. Mit diesem Vektor k^ gilt 

p(§L w + h w -l£ T ) <1 (25) 

Bezeichne nunmehr A* einen Eigenwert von — • $ + — ■ h ■ k*T 

und \ (wie oben) einen Eigenwert von & w := ^ — h-k^f sowie E_ die Einheitsmatrix, so folgt aus 
= det ( X* ■ E - ( — ■ $ + — ■ h ■ k*J ) ) = det ( w ■ X* ■ E - ($ + h ■ k* J 



s w w 
endlich Xi — w ■ X* . Mit diesem fc^, gilt: 

P 



< w (26) 



Das so erweiterte Iterationsverfahren wird mit den Rechnungen des nachfolgenden Unterpunkts 
vorgestellt. Es wird wieder das unter 1.3 vorgestellte Rechenbeispiel aufgegriffen; diesmal aber mit 
9 bzw. 12 Nachkommastellen. 
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3.2 Rechenbeispiele zum erweiterten Iterationsverfahren 



Mit 



59 


-63 


-2 


29 


42 


51 


3G 


31 


-67 




A = 29 42 51 und 



berechnen sich die $ - Matrix und der h - Vektor wie in 1.3. 
Mit w = jq und $„, = i • $, h w = ± ■ h lauten 

0,000.000.000 10,677.966.102 0, 338.983.05l\ / -12,372.881.355.932^ 

$ w = | -6,904.761.905 0,000.000.000 -12,142.857.143 und h w = 63,333.333.333.333 

5,373.134.328 4,626.865.672 0,000.000.000 / \ 15,373.134.328.358, 

Unter Verwendung von S_, c und "r" wie in 2.2 erhält man für N = 10 und der Einheitsmatrix für 
P w mit Gleichung (20): 



'48,959.414.201 29,945.482.078 26, 524.385. 209^ 
23,945.482.078 150,581.701.310 12,904.368.225 
,26,524.385.209 12,904.368.225 36,708.746.563; 



/ 6.195, 608.211.100 30.591,456.824.000 -774,625.645.620 

P x ss 30.591,456.824.000 181.987,753.000.000 -7.760,651.406.800 
V -774,625.645.620 -7.760,651.406.800 622,611.391.310 

Für die Vektoren fc^ ergibt das mit Gleichung (19): 

kj « (-0, 080.594.648.281 - 0, 013.857.607.126 - 0, 175.695.545.308) 



fco « (-0, 117.736.653.973 0, 026.270.720.999 - 0, 199.437.647.693) 



Mit k* w = k Q folgt sodann mit: 



für 



*T 



1™ := £ - ä ■ k w 



und für 



A ■ x r 



und den Startvektor x£ = (1,1,1) die Vektoren: 

2,971.843.010.215\ / 1, 000.001. 061. 532\ / 0, 000.048.722 N 

m + l = 5: 5,943.680.457.527 x 5 = 2,000.000.251.191 r 5 = -0,000.008.825 

5,915.516.643.666/ \ 2,999.999.016.477/ \ 0,000.111.898, 



1 = 10 



2, 971.840.542. 360\ 
5,943.681.084.719 x w 
8,915.521.627.079/ 



1, 000.000.000. 000^ 

2,000.000.000.000 

3,000.000.000.000; 



LW 



'-1,421.085 10~ 14N 
0,000.000 io~ 14 
-2,842.171 10~ 14 
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Die Mächtigkeit des erweiterten Iterationsverfahrens zeigt sich erst bei Erniedrigung des Parame- 
ters "w" ! 



Seien w = und die anderen Parameter unverändert, dann erhalten wir nach N = 5 Iterations- 
schritten (beginnend mit P 5 = Einheitsmatrix): 

4, 877.396.786 10 7 2, 975.175.473 10 8 -1, 334.071.540 10 7N 
P x « | 2,975.175.473 10 8 1, 818.776.612 10 9 -8.187.822.688 10 7 
-1, 334.071.540 10 7 -8.187.822.688 10 7 3, 712.634.543 10 6 



fco « (-0, 117.898.057.028 0, 025.378.263.771 - 0, 199.405.203.705) 

und nach M = m+1 = 5 Schritten für die Fixpunktiteration 

2,988.260.385.409\ / 1, 000.000.000.000\ / 4,908 10~ ir 

M = 5: x 5 =\ 5,976.520.770.812 x 5 = 2,000.000.000.000 r 5 = -8,868 10~ 12 



8,964.781.156.214/ \ 3,000.000.000.000/ \ 1,112 10 
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Seien jetzt noch w = 10 o 00 , N = 2 und M = m+1 = 4, dann erhalten wir (allen anderen Parameter 
unverändert): 

/ 2,061.646.012 10 14 1, 260.667.803 10 15 -5,678.165.532 10 13N 
P ± « 1,260.667.803 10 15 7, 708.847.447 10 15 -3, 472.169.309 10 14 
\-5, 678.165.532 10 13 -3, 472.169.309 10 14 1, 563.938.679 10 13 

fco ks (-0, 117.899.526.977 0, 025.369.284.637 - 0, 199.404.800.048) 



2,988.423.068.799\ / 1, 000.000.000.000N 

x A =\ 5,976.846.137.597 = 2,000.000.000.000 r 4 = 

8,965.269.206.396/ \ 3,000.000.000.000/ 

3.3 Algorithmische Formulierung des erweiterten Iterationsverfahrens 

Vorbereitend für die algorithmische Formulierung werden die folgenden Definitionen eingeführt: 
Definition 3.1 Die Funktion Rice 
Rice : R nxn — ► R" xn 



7, 958 


io- 


-13 


1,705 


io- 


-13 


1,691 


10" 


-12 



P^Ricc(P) = Q + ^-U-P-(k w .hl)-P 



1 



Definition 3.2 Die Funktion F_ 

F : R™ — > R™ 



Nun zum Algorithmus. 
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1. Zum linearen Glcichungssystcm A ■ x — b : $ und h berechnen. 

2. < w < 1 , N und M wählen. 

3. $„„ := — ■ $ und h,„ :— — ■ h berechnen. 

4. Mit Startwert P_ N = E - Einheitsmatrix in N - Iterationsschritten mittels Funktion "Rice" 
durch 

£„_! =Ricc{P n ) 

den Näherungswert P_ x berechnen. 



5. Den Vektor k berechnen: 



6. Mit beliebigem Startwert x in M - Iterationsschritten (M > 1) mittels der Funktion "F_ 
durch 

äm+i =E(£.m)+ä 

den Näherungswert x M berechnen. 

7. Den Näherungswert x M gemäß Vorschrift: 

_( 1 \ - 

und 

Ha/ = 4 • £m — 

berechnen. 
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